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Abstract 

It is shown that the problem of calculating spin-spin correlation functions, 
in the dimers RVB states, on a possibly diluted 2D square lattice, can be 
formulated in terms of a transfer matrix. The transfer matrix is used for 
exact numerical calculations of spin-spin correlation functions on ladders up 
to four units wide. 
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In the last few years, the class of Resonating Valence Bonds (RVB) states have drawn 
much attention in connection to the long lasting pursuit of the anti-ferromagnetic 2D Heisen- 
berg model ground state [0,0. Especially important is the case of spin half, which is con- 
sidered to be relevant to high Tc superconductors. In this context, simulations had been 
performed to estimate the spin-spin correlation functions in these states || ; and particularly 
the expectation values of the spin half Heisenberg Hamiltonian were calculated for varia- 
tional considerations. In the two legged ladder lattice the short range RVB states, which are 
in the focus of this article, were found to have low energies @]. Despite this interest these 
simulations had not been backed by exact results. This article provides a way to extract such 
results, in the case of the 2D dimers RVB states on a (possibly) diluted lattice, by means of 
reducing the problem of calculating the spin-spin correlation functions, to the considerably 
simpler problem of treating a transfer matrix defined on a ID lattice. 

On a general lattice the class of RVB states is defined by 

i*> = £ n ua,3)(uij-yii) c 1 ) 

{n}(ij') e n 

where {11} are all the possible divisions of the lattice into pairs of sites - in which any 

site is a member of one and only one pair, and U (i, j) is nonnegative. For the short-range 
RVB states U has a cutoff. We will restrict | ty) to be a dimers RVB state on a square 
lattice, hence the function U is given by 

x if i — j = x 

U = \ y if i-j = y (2) 

otherwise 

where x and y are unit vectors. In this case we can replace {11} by the dimers coverings of 
the lattice, denoted by {A}. 

The spin-spin correlation function in the dimers RVB states is defined by the expectation 
value |,|] 



where i and j are two sites on the lattice. Two quantities are evaluated for the calculation 
of this expectation value. The norm : 



{&lAr} 

= y MK J2 2**71"" (4) 
{AlAr} 

where A^ and A# are any two dimers coverings of the lattice which are placed on each other 

to create {A^, Ar} which is the ensemble of all the loops configurations given by overlaps 

of dimers coverings, this ensemble on a - 2 x 2 lattice is depicted in Fig. |l|; n x and n y are 

the numbers of horizontally and vertically placed dimers respectively (with n x + n y = N, 

where N is the number of sites in the lattice); A/a is the number of loops in this overlap; 

and 77 = -. And 
' y 

c ij = ±l^\s i -s j m 

= y MK £ 2 M Xr] n x (5) 

where the sign is '+' when the two sites are on the same sub-lattice and '— ' otherwise, and 
{A^, Af>}^j is the ensemble of all the loops configurations in which the sites i and j are on 
the same loop - X^j (this loop is particular to each of the loops configurations). See Fig. [I] 
for further explanation. 

As an introduction, let us recall the construction of the transfer matrix, formulated to 
solve the dimers problem |J, that is, to calculate the partition function 

z = Y J y ny x nx . (6) 

{A} 

In an M rows by K columns square lattice the transfer matrix - V is an operator defined 
in a 2 K dimensional Hilbert space of the Ising states of K spins half placed in the sites of 
a K sites row, each in a site. In a certain dimers configuration, each row in the lattice is 
represented by an Ising state. This Ising state denotes the vertical dimers which are placed 
on the bonds between this row and the row above it. An up spin in a site signifies the 



presence of a vertical dimer on the bond between this site and the site above it (an up- 
dimer), while a down spin signifies the absence of an up-dimer on this bond. We will use the 
Fock representation where the spins are represented by Bosons (or Fermions) of spin half 
and 

\si, ..., SK )=n(4 Si 10)*) (7) 

i 

where s* =f or j. 

Given an up-dimers configuration of a certain row, starting with the first row, V will 
generate all the permitted dimers configurations of the next row, each marked by its up- 
dimers. V is copmposed of two operators : 

V = V 2 V X . 

The first is 

vi = n <^ (8) 

1=1 

where 

This operator reverses all the spins for the next row. After this operation a vertical or 
horizontal dimer can be placed on any up-spin site. A vertical dimer is placed by leaving 
the spin up in the site, and a horizontal dimer, between any two adjacent up-spin sites % and 
+ is placed by the operator a^a^ +1 . The operator which places m horizontal dimers in 
arbitrary locations along the row is (X^=i a i a i+i) ■ An arbitrary number of horizontal 
dimers, each accompanied by the weight rj, are placed by the operator 

K 



V 2 = exp ^ Yl a i a m) 



(9) 



where in the case of periodic boundary conditions (which are henceforth assumed) <J K+1 
erf . Hence we get (note that in Eq. (||) n x + n y = MK/2) 



Z = y l 2 MK TrV M . (10) 
Let us now consider the square of Eq. (Q) 

Z 2 = y ny x nx . (11) 

{A l ,Ar} 

This is just summing over loops configurations in whice each site is connected by dimers 
of two flavours - L and R, which are independent of each other. To formulate a transfer 
matrix, we denote each of the flavoured dimers in any site, by a flavoured spin. These spins 
are propagating along the columns by the operator 

v = V L V R 

where Vl and Vr are defined in a Hilbert space which is the tensor product of a L and R 
Ising spaces, and Vl (Vr) effects only the L (R) spins in each of the product states. 

After the prepatatory examples above we are ready to tackle the transfer matrix repre- 
santation of Eq. (^), with the additional term of 2^ x for each loops configuration. Imagine 
that there are two colours of dimers - red (r) and green (g). A dimer from each of the colours 
is distincted by an additional flavour index - a £ {L, R}. We will define the ensemble {0} 
to include all the possible coloured-dimers configurations in which each site in the lattice is 
connected with two dimers, of identical colours, marked by distinct flavour index. In {0} 
there are exactly 2 Jsfx different configurations for any loops configuration in Eq. ([TTJ) , since 
each of the loops in this overlap can appear in two forms - all red, or all green; and we can 
write 

V = y MK Y,V nx - (12) 

{6} 

We denote the coloured dimers in each site by coloured spins. Accordingly we expand 
further our Hilbert space to be the 8 K dimensional space of a row of K sites; on each site 
- i two spin half Bosons, of the flavours L and R, both labeled by the same colour index - 
Q £ {r, g}. Each of the states in this space is specified by 
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I S (1,L); S (l,-R); c l> • • • > S (K,L), S (K,R), C K ) = 
K 

II ( C i,Ci,L,a (i , L) C i,Ci,R,s {i>R) ( 13 ) 

i 

Flipping spins must now include the two colour options so 

K 

Vl = II (<r,L<r,H + ( 14 ) 
i=l 

where for example 

Horizontal dimers are placed by lowering spins of the same flavour and colour by four 
operators of the form 

V^ a = exp ^ c ,«^ + i, C) „j (15) 

where a = L,R and c — r,g, so now 

v 2 = n v? a - (i6) 

a = L, R 
c = r } g 

But this is not all. The operator defined in Eq. (fPf ) preserves the colour of a column, 
hence imposing the mono-colour propagation of a loop with a vertical dimer. But if a site is 
not connected by a vertical dimer from 'below', its pair of dimers should be of either colours, 
as their are two colour options for the loop passing through this site. Since the absence of 
vertical dimers is denoted by two down spins, the additional colour flipping options are 
provided by the operator 

K 

Vo = J] (l + C \,r,L,l C Xr,R,l C i,9,L,lCi, g ,R,l 

i=l 

+ C i,g,L,i C i,g,R,i Ci ^ L ri C i,r,R,i) • ( i7 ) 

Combining equations (@),(0), and (fPTD, the transfer operator is 
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V = V 2 ViV (18) 

and 

V = y MK TrV M . (19) 

On narrow lattices, such as the two legged ladder, it is possible to extend this treatment to 
longer range, like the first and second, or the first and third nearest neighbours, RVB states 
0. This extension, although straightforward, requires a more complicated transfer matrix, 
which is usualy larger in its dimesion (for a given lattice width). 

The transfer matrix formulation will be concluded in adding, that a disconnected site (a 
static hole) i, can be placed by positioning the operator 

Vf = - <tT ct7 (20) 

* y l,r,L l,r,R y ' 

to the left of Vi, when propagating in to i's row. This operator colours the disconnected site 
in red. A class of variational states, that include annealed holes, which obey Fermi statistics, 
will be treated elsewhere, using a similar technique. 

In order to calculate C^j defined in Eq. (pj), we will calculate the quantity y^j, defined 

by 

yij = y MK E 2A V* (21) 

where {^}ij is the ensemble of all the loops configurations in which the sites i and j are not 
on the same loop. Some of these configurations are easily counted, by fixing the two pairs of 
dimers, placed on the two sites, to be of distinct colours, or by calculating the colour-colour 
correlation function 

y A= MKrj, r (y{M-p) r yp g \ (22) 

where the two sites are p rows apart, and 

n i, a = n L, T + n h,l 



(the flavour index - a is arbitrary). In Eq. ( p2|) we calculate the contributions to y^j 
from configurations in which A^ is a green loop and A^- is a red loop. Each of the other 
configurations in which Aj 7^ A^- (that is, loops configurations which do not contribute to 
y^jj but must be counted in y^j) has its duplicate in one and only one of the configurations 
which do contribute to 3^J\ since each of the two loops may apear in either of the two 
colours. Because there are four colours possibilities for the two loops, we conclude that 

3>*i = 4 • y%. (23) 

Finally, since the sum for the norm in Eq. @ contains all the loops configurations, those 
in which the two sites are on the same loop, and those in which they are not, we get 

C ij= V - y ij- ( 24 ) 

The spin-spin correlation functions were numericaly calculated, in the isotropic dimers 
RVB state (y = rj = 1), on non-diluted ladders of the sizes 2,3,4 x 40, with vertical (along 
the ladders legs) periodic boundary conditions, and horizontal periodic boundary conditions 
in the case of 4 x 40 lattice. The results are summerized in table |IJ. 

In conclusion, using the transfer matrix method it was shown, that it is possible to 
reduce the complexity of calculating spin-spin correlation functions in dimers RVB states 
on possibly diluted 2D lattices, to the complexity of a ID quantum many-body problem. 
This simplification permited exact calculation of correlation functions on narrow lattices. 
We suggest that these kind of calculations are possible when longer range RVB states are 
considered. 

The author thanks Assa Auerbach for essential guidance. This work was supported by a 
grant from the Israely Academy of Science, and the Foundation for Promotion of Research 
at Technion. 



S 



REFERENCES 

* Email: havilio@pharaoh.technion.ac.il 

[1] P.W. Anderson, Science 235, 1196 (1987). 

[2] For a general review see : Interacting Electrons and Quantum Magnetism, Assa Auerbach 
( Springer- Verlag, New York, 1994), ch. 8. 

[3] S. Liang, B. Doucot, and P. W. Anderson, Phys. Rev. Lett. 61, 365 (1988). 

[4] S. R. White et al, Phys. Rev. Lett. 73, 886 (1994). 

[5] M. Kohmoto and Y. Shapir, Phys. Rev. B 37, 9439 (1988). 

[6] Elliott H. Lieb, J. Math. Phys. 8, 2339 (1967). 

[7] Moshe Havilio unpublished. 

[8] This result, for the correlation length on the two-legged ladder, disagree with the result 
in Ref. 4. 

[9] In the case of 3 units width lattice, the correlation function is slightly wavy, hence a 
slightly different result may apear, if the correlation length is to be calculeted by the 
values of the correlation function in different distances. 



9 



TABLES 



width 


2 


3 


4 


(a) 


0.286139 


0.328752 


0.246393 


(b) 




0.277263 




(c) 


0.539779 


0.344843 


0.395096 


\E\ 


0.556029 


0.656432 


0.641489 


(a) 


0.696652 @ 


0.735272 @ 


0.723966 


(b) 




0.780763 





TABLE I. Results of the one lattice unit spin-spin correlation functions, energies, and corre- 
lation lengthes, for the isotropic dimers RVB state on non-diluted ladders of the width 2, 3, and 
4 lattice units. The results in the second row are for the absolute values of the one lattice unit 
correlations: (a) in the vertical case along one of the side legs of the ladders; (b) in the vertical 
case along the middle leg of a 3 x 40 ladder; and (c) in the horizontal case. The results in the third 
row are for the expactation values of the isotropic Heisenberg Hamiltonian for each of the lattices. 
The expactation values are in units of J /site. The results in the fourth row are for the correlation 
lengthes along the (a) side legs of the ladders, and (b) middle leg of a 3 x 40 ladder. The results 
for the correlation lengthes are according to the values of the correlation function in distances of 
five and six lattice units along the lattice. 
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FIG. 1. The four possible loops configurations, made by overlaps of dimers coverings, on a - 
2x2 lattice. The dimers configurations from the L side are placed on the dimers configurations 
from the R side. In the sum of Eq. (5), for the two sites marked by i and j, the first (from left to 
right) three loops configurations have to be included. 
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